Slow 4 He Quenches Produce Fuzzy, Transient Vortices 
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We examine the Zurek scenario for the production of vortices in quenches of liquid 4 He in the light 
of recent experiments. Extending our previous results to later times, we argue that short wavelength 
thermal fluctuations make vortices poorly defined until after the transition has occurred. Further, 
if and when vortices appear, it is plausible that that they will decay faster than anticipated from 
turbulence experiments, irrespective of quench rates. 
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Several years ago Kibble initiated a programme jlj] for 
examining the consequences of phase transitions in the 
early universe for large scale structure formation. In 
particular, cosmic strings produced at these transitions 
were seen to have the potential for creating the observed 
fluctuations in the microwave background. Using simple 
causal arguments [|| he obtained bounds on the density 
of strings produced in such transitions. Unfortunately 
such consequences can only be inferred indirectly and 
with ambiguity. On the other hand, the production and 
measurement of vortices in condensed matter (the coun- 
terpart of cosmic strings in OFT) is a viable experiment. 
In a series of papers Zurek || made predictions for the 
density of vortices produced during quenches of liquid 
4 He and 3 He into their superfluid states, using similar 
causal arguments. 

Zurek's original assumption was that, as we approach 
the (continuous) transition, the correlation length £(t) 
of the field initially matches the equilibrium correlation 
length £ eq (t) of the field. However, once is growing 
at the speed of (second) soundF], at some time t before 
the transition begins (as characterised by the divergence 
of £eq(i)) it freezes. Only after the transition has begun 
does the field unfreeze, at a comparable time interval i. 

With superfluid He in mind we assume a time- 
dependent Landau- Ginzburg free energy of the form 



F(t) = 



d 3 3 
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In F(t), (f> = (tfii + i<p 2 )/V2 is the complex order- 
parameter field, whose magnitude determines the su- 
perfluid density. In equilibrium at temperature T, 
in a mean field approximation, the chemical potential 
a(T/T c ) takes the form a(T/T c ) = a e(T/T c ), where 
e = (T/T c — 1) measures the critical temperature T c rel- 
ative to T. At a quench the change in T/T c with time 
leads most simply to e of the form 




*The corresponding argument in QFT replaces the speed of 
sound by the speed of light. 



e(t) = e 0(t) 

T Q 

for — oo < t < tq(1 + eo), after which e(t) — — 1. eo 
measures the original relative temperature and tq defines 
the quench time. The quench begins at time t — and 
the transition from the normal to the superfluid phase 
begins at time t = to = £ot~q- 

With £q = h 2 /2mao and To = h/ctQ setting the funda- 
mental distance and time scales, mean-field critical be- 
haviour fixes H i = ^/tqtq. The frozen correlation length 
of the field £ = eo(T Q /T ) 1 /4 ~ £ eq (t - t) « £ eg (i + t). 
The ordering of the field on large scales will be marred 
by topological defects, in this case vortices. If the ini- 
tial vortex separation is £def , then their density, Udef , is 
n def — 0(l/£jL/)- Zurek initially made the second as- 
sumption (matched by Kibble in QFT) that £def = 0(£) 
whereby 



ndef 



where / = O(l). Since £o also measures cold vortex 
thickness, tq 3> tq corresponds to a measurably large 
number of widely separated vortices. 

This is a remarkable prediction in that it says that, 
in the first instance, the dimensionful interactions /3|^| 4 
introduce no new scale. It has been the motive for several 
experiments. One p| at Lancaster on vortex production 
in He shows consistency with the prediction, whereas a 
second experiment 0j on 4 He finds no vortices. On the 
other hand, vortex production in quenches of 3 He into 
its B-phase in two separate experiments at Grenoble || 
and Helsinki |?J show quantitative agreement with ndef 
above. This is even though it is now argued || that the 
Helsinki experiment should not show agreement because 
of the geometry of the heating event. 

In an earlier letter ||] we argued that, RG improvement 
apart, the simple prediction above for ndef is unreliable 
on two counts. Firstly, and crucially, vortex separation 
cannot be related simply to field or phase correlation 
length. It is much more sensible to relate vortex den- 
sity to the density of line zeroes, n zero , of the field, since 
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a classical vortex is a tube of false vacuum, with a line 
zero at its core. For Gaussian field fluctuations 

(0„(r,i)&(O,t)) = 5 ab G{r,t), 



the line-zero density n zero (t) can be shown ]f 0|Jlf[ | to have 
the form (where primes denote differentiation with re- 
spect to r) 



ro(t) 



-1 G//(0,t) 
2?r G(0,t) ' 



It depends manifestly on the short- distance behaviour of 
G(r, t) , rather than the correlation length £ that follows 
from large-distance behaviour. Secondly, if a line zero is 
to be counted as a vortex, at the very least it should have 
fractal dimension close to unity on a scale comparable 
to the classical vortex thickness £oj even though it may 
behave like a random walk on larger scales. Otherwise 
the density would depend on the scale at which we look. 
This is equally true for cosmic strings [H . 

In P] we analysed the model until the time that the 
transition was effected. In this letter we shall now ar- 
gue that, even when the transition is complete, for slow 
quenches the short wavelength thermal fluctuations are 
likely to persist for some time at a level that prevents an 
approximately scale-invariant definition of vortex den- 
sity. Further, after that time it is plausible that such 
vortices as do exist decay more rapidly than is assumed 
in the analysis of ||. 

Motivated by Zurek's later numerical 0,^5| simula- 
tions, we adopt the time-dependent Landau-Ginzburg 
(TDLG) equation for F, 



1 d(f> a 
T dt 



SF 



where r\ a is Gaussian thermal noise, satisfying 

(?? a (x, t) Vb (y', f )) = 2TT(t)6 ab S(x - y)S(t - t'). 

In space, time and temperature units in which £o — 
To = fc_B = 1, the TDLG equation becomes 

a (x, t) = -[-V 2 + e(t) + 0|0(x, i)| 2 ]av(x, t) + 77a(x, t), 

where (3 is the rescaled coupling and fj the rescaled noise. 
We adopt the approximation of preserving Gaussian fluc- 
tuations by linearising the self-interaction as 



^a(x,t) 



e//W]^a(x,t) + 77 (x,t), 



where e e // contains a (self-consistent) term 0(f3(\cf>\ 2 )). 
Additive renormalisation is necessary, so that e e ff w e, 
as given earlier, for t < to- 

Self-consistent linearisation is the standard approxi- 
mation in non-equilibrium QFT fTq| , but is not strictly 
necessary here, since numerical simulations that identify 
line zeroes of the field can be made that use the full self- 
interaction |l5| ]. However, to date none address the ques- 
tions we are posing here exactly, and until then there is 



virtue in analytic approximations provided they are not 
taken too seriously. 

The solution for G(r, t) from the linear equation for <p a 
is straightforward, most usefully in the Schwinger proper- 
time representation, as G(r, t) = 



drTOt-T/i)'. ' 

V 47TT 



3/2 



-r 2 /4r - Jo ds e =//(*- s /2) 



where T is the rescaled temperature. 

For constant e e //, as happens at early times, we re- 
cover the usual Yukawa correlator. However, a saddle- 
point calculation at time t n gives £(fo) = 0(£), confirm- 
ing Zurek's causal result for the freezing of the field cor- 
relation length. 

A direct substitution of G(r, t) into n zero is impossi- 
ble because of UV singularities at small r. However, 
a starting-point for counting vortices in superfluids is 
to count line zeroes of an appropriately coarse-grained 
field, in which structure on a scale smaller than £oj 
the classical vortex size, is not present. This is, in- 
deed, the basis of Zurek's and other numerical simula- 
tions. For the moment, we put in a momentum cutoff 
k^ 1 > I = Z£o = O(£o) by hand, which corresponds to 
damping the singularity in G(r, t) at r = as 

r°° dr T(t - t/2) 2 /4t £ ds eeff(t - s /2) 
° U ' j Jo [Mr + H 3/2 

making G;(0,i) finite. We stress that, for t as to, the 
correlation length £ remains 0(£)j independent of I. 

Assuming a single zero of e e / / (t) at t = to, at r = the 
exponential in the integrand peaks at t = f = 2(i — to). 
Expanding about f to quadratic order gives 



Gi(0,t)*sT c e J *o ' in 



[47r(r + r 2 )] 3 / 2 ' 



The effect of the back-reaction is to stop the growth of 



Gi(0,*) - G,(0,to) 



)o at its symmetry- 



broken value in our dimensionless units. We see that 
this cessation of growth requires lim^^oo e e ff (u) = 0. 
That is, we must choose 

e e //(t) = e(t)+ / 9(G l (0,t)-G,(to s 0)), 

thereby preserving Goldstone's theorem. 

Beyond that, what is remarkable in this approxima- 
tion is that the density of line zeroes uses no property 
of the self-mass contribution to e e f / (t) , self-consistent or 
otherwise. In n zero all prefactors cancef], to give 



n>zero [t) 



1 r°° dr ( T _2(t-t )) 2 /4P 

1 JO ( r+ ;2)5/2 c 



47T f°° iL e -(r-2(i- 

JO ( T +P)3/2 e 



t ))V4t 2 



^Our ignoring prefactors in |j] was fortuitous, leaving our 
conclusions obtained there unaffected. 
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on using the definition tq = P in natural units. At t = to 
both numerator and denominator are dominated by the 
short wavelength fluctuations at small r. Even though 
the field is correlated over a distance £ 3> I the density of 
line zeroes n zero — 0(l~ 2 ) depends entirely on the scale 
at which we look. However, as time passes the peak of 
the exponential grows and n zero becomes increasingly in- 
sensitive to I. How much time we have depends on the 
magnitude of /?, since once G(0, t) has reached this value 
it stops growing. Since G(0, t) = 0(exp(((f — to)/t) 2 ) at 
early times the backreaction is implemented extremely 
rapidly. We can estimate the time t* at which this hap- 
pens by substituting e(it) for e e //(u) in the expression for 
Gi(Q,t) above. 

For t > t* the equation for n zero (t) is not so simple 
since the estimate above, based on a single isolated zero 
of e e t f (t) , breaks down because of the approximate van- 
ishing of e e f f (t) for t > t*. A more careful analysis shows 
that Gi(0,t) can be written as 



Gz(0,t) 



drT{t-T/2) 
[4tt(t + P)} 3 / 2 



G(r,i), 



where G(r, t) has the same peak as before at r = 2{t— to), 
in the vicinity of which 

G{r,t) = e 2 ft du{ee " (u)l e-i-r-W-to)) 2 /*? ^ 

but G(r,i) = 1 for t < 2(t - t*). Thus, for t q > r , 
Gi(0,t) can be approximately separated as 



G l (0,t)^GY v (t) + G lu {t) 



IRi 



Gf v (t)=f(t) / dr/[47T(r + P)f 2 



where 



describes the scale- dependent short wavelength thermal 
noise, proportional to temperature, and 



G 1H {t) 



(8tt(* - t )) 3 / 2 



drG{T, t) 



describes the scale-independent, temperature indepen- 
dent, long wavelength fluctuations. A similar decompo- 
sition G// ; (0,t) ^ G//Y v {t) + Gn IR {t) can be performed 
in an obvious way, with G// IR (t)/G IR (t) =0(t~ 1 ). 

Firstly suppose that, for t > t* , G IR {t) > Gf F (t) 
and Gn IR (t) > Gf>Y v (t), as would be the case for a 
temperature quench T(t) — ► 0. Then, with little thermal 
noise, we have widely separated line zeroes, with density 
n de f(t) « -G// IR (t)/2TrG IR (t). With dn zero /dl small in 
comparison to n zero /l at I = £o we identify such essen- 
tially non-fractal line-zeroes with vortices, and n zero with 
ndef- Of course, we require non-Gaussianity to create 
true classical energy profiles. Nonetheless, the Halperin- 
Mazenko result may be well approximated even when the 



fluctuations are no longer Gaussian Wm. This is sup- 
ported by the observation that, once the line zeroes have 
straightened on small scales at t > t*, the Gaussian field 
energy, largely in field gradients, is 

Pk (//^(W 2 ) =-VG"(0,t), 
where V is the spatial volume. This matches the energy 
E » Vn def (t)(2TrG(0,t)) = -VG"(0,t) 

possessed by a network of classical global strings with 
density n zero , in the same approximation of cutting off 
their logarithmic tails. 

For times t > t* , but not too late, 



n-def(t) 



t 




the solution to Vinen's equation fil 



dn 



def 



Ot 



-X2—n def , 
m 1 



where \2 — 47r is a universal constant. 

This decay law is assumed in the analysis of the Lan- 
caster experiments, in which the density of vortices is in- 
ferred from the intensity of the signal of scattered second 
sound. RG improvement leaves 3 He unchanged, but for 
4 He it redefines X2 to x 2 (l -T/TJ- 1 / 3 > X 2- The prob- 
lem is that this makes an already large \2 even larger. 
In the attenuation of second sound the signal to noise 
ratio is approximately 0(l/x2*)- The empirical value of 
X2 used in the Lancaster experiments is not taken from 
quenches, but turbulent flow experiments. It is suggested 
|| that X2 ~ 0.005, a good three orders of magnitude 
smaller than our prediction above. Although the TDLG 
theory is not very reliable for 4 He, if our estimate is 
sensible it does imply that vortices produced in a tem- 
perature quench decay much faster than those produced 
in turbulence. 3 He experiments provide no check. 

In reality, the situation for the Lancaster A He exper- 
iments is even more complex, since they are pressure 
quenches for which the temperature T is almost con- 
stant at T ~ T c . Unlike temperature quenches Jl5| , ^6[ , 
thermal fluctuations here remain at full strength. Even 
for 3 He, T/T c never gets very small, and henceforth we 
take T = T c in G;(0,t) above. The necessary time- 
independence of G IR (t) for t > t* is achieved by taking 
e e ff(u) = 0(u~ l ). In consequence, as t increases beyond 
t* the relative magnitude of the UV and IR contribu- 
tions to G;(0,t) remains approximately constant. Fur- 
ther, since for t = t* , 



a a/ t * du\e ef f{u)\ (At) 2 /* 2 



1. 



this ratio is the ratio at t = t* . 

Nonetheless, as long as the UV fluctuations arc in- 
significant at t = t* the density of line zeroes will 
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remain largely independent of scale. This follows if 
Gtt IR (t*) > GtlY v [t*), since G//;(0,i) becomes scale- 
independent later than G/(0, i). In M we showed that 
this is true provided 



(tq/t )(1-T g /T c ) <Ctt\ 

where C = O(l) and Tq is the Ginzburg temperature. 
With tq/t = O(10 3 ) and (1 - T G /T C ) = O(10~ 12 ) this 
inequality is well satisfied for a linearised TDLG theory 
for 3 He derived^ from the full TDGL theory gjjj, but 
there is no way that it can be satisfied for 4 iJe, when 
subjected to a slow mechanical quench, as in the Lan- 
caster experiment, for which tq/tq = O(10 10 ), since the 
Ginzburg regime is so large that (1 — Tq/T c ) = O(l). 
Effectively, the 4 He quench is nineteen orders of mag- 
nitude slower than its 3 He counterpart. It is satisfying 
to find the Ginzburg regime, originally pp9|] but erro- 
neously thought to determine the initial defect density, 
but then discarded |^^|, reappear as an essential ingre- 
dient in undermining universality. 

Thus, for 3 He, with negligable UV contributions, we 
estimate the primordial density of vortices as 



the earlier problem of the large X2 = 0(f 2 ). In the ab- 
sence of any mechanism to reduce its value drastically, 
this would make it impossible to see vortices. As a sep- 
arate observation, we note that the large value of f 2 in 
the prefactor of n^e/ is, in itself, almost enough to make 
it impossible to see vortices in 4 He experiments, should 
they be present. This will be pursued elsewhere. 

In summary, this work suggests that for slow pressure 
quenches in 4 He, we see no well-defined vortices at early 
times because of thermal fluctuations, and it is plausible 
that, if we do see them at later times, there are less than 
we would have expected because of their rapid decay and 
their initial low density. The situation is different for 
3 He. That the density of such vortices as appear may 
agree with the Zurek prediction is essentially a conse- 
quence of dimensional analysis, given that the main ef- 
fects of the self-interaction have a tendency to cancel in 
the counting of vortices, without introducing new scales. 

This work is the result of a network supported by the 
European Science Foundation. I thank Tom Kibble, Gr- 
isha Volovik, Nuno Antunes and Eleftheria Kavoussanaki 
for helpful conversations and analysis. 



n def {t*) 



8vr(r - t ) e 2 




in accord with the original prediction of Zurek. Because 
of the rapid growth of G_(0, t), (t*-t )/i=p> 1 = 0(1). 
With p behaving as (In ft) 1 / 2 there is very little variation. 
For 3 He quenches p sa 5 (and for 3 He quenches p ~ 3). 
We note that the factor^ of f 2 = 8irp gives a value of 
/ = 0(10), in agreement with the empirical results of || 
and the numerical results of 0,|h|q 

However, if the inequality is badly violated, as hap- 
pens with 4 ife for slow pressure (but not temperature) 
quenches, then the density of zeroes ridef = 0{l~ 2 ) after 
t* depends exactly on the scale I at which we look and 
they are not candidates for vortices. Since the whole of 
the quench takes place within the Ginzburg regime this is 
not implausible. However, it is possible that, even though 
the thermal noise never switches off, there is no more 
than a postponement of vortex production, since our ap- 
proximations must break down at some stage. The best 
outcome is to assume that the effect of the thermal fluc- 
tuations on fractal behaviour is diminished, only leading 
to a delay in the time at which vortices finally appear. 
Even if we suppose that ndef above is a starting point 
for calculating the density at later times, albeit with a 
different to, thereby preserving Vinen's law, we then have 



^Ignoring the position-dependent temperature of 
^An errant factor of 3 appeared in the result of [ 
"The temperature quench of the latter is somewhat differ- 
ent from that considered here, but should still give the same 
results in this case 
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